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Abstract 

Point-like Liouville integrable dynamical defects are introduced in the context of the 
Landau-Lifshitz and Principal Chiral (Faddeev-Reshetikhin) models. Based primarily on the 
underlying quadratic algebra we identify the first local integrals of motion, the associated 
Lax pairs as well as the relevant sewing conditions around the defect point. The involution 
of the integrals of motion is shown taking into account the sewing conditions. 
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1 Introduction 

A systematic algebraic formulation for describing Liouville integrable point-like defects was 
recently introduced [1, 2]. The description was primarily based on the underlying quadratic 
algebra satisfied by the bulk monodromy matrices, which describe the left and right theories, 
as well as by the defect L- matrix. In fact, this is the main necessary requirement so that 
Liouville integrability may be by construction guaranteed. The main steps of the proposed 
algebraic methodology are provided in the subsequent sections, however for a more detailed 
description of the process we refer the interested reader to [1, 2]. 

The issue of integrable defects has been quite intriguing, and a considerable number of 
studies have been devoted to this particular problem both at the classical and quantum 
level [3]-[18]. The problem becomes even more intricate at the level of classical integrable 
field theories, where usually the defect is introduced as a discontinuity (jump) together with 
suitable sewing conditions [6, 8]. In the present scheme the sewing conditions naturally 
emerge as continuity conditions on the time components of the Lax pair around the defect 
point. It was also rigorously proven in [1] that the sewing conditions are compatible with 
the hierarchy of the Hamiltonians, a fact that ensures that the proposed formulation is well 
defined and consistent. 
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In the present work, we consider dynamical point-like defects in the context of sigma 
models, such as the Landau-Lifshitz (L-L) model and a variation of the familiar principal 
chiral model (PCM), the so called Faddeev-Reshetikhin (F-R) [19] model, in such a way that 
integrability is preserved. It is worth noting that this kind of systems, in addition to their 
own physical and mathematical value have also attracted considerable interest lately due 
to the fact that they typically arise within the AdS/CFT context (see e.g. [20, 21, 22] and 
references therein). Thus, it would be of great significance to investigate possible relevant 
physical implications in this particular frame. Recall also that the typical PCM possesses a 
non-ultra local algebra rendering its quantization a very intricate task, whereas its variation, 
the F-R model is associated to a more familiar ultra-local algebra, much easier to deal with. 
In any case, both PCM and F-R models share the same Lax pair, and consequently the 
same equations of motion, hence they are physically quite similar. Therefore, the findings 
presented for the F-R model are naturally relevant for the conventional PCM model. 

Based on the formulation of [1, 2] we introduce the defect matrix, and the associated 
modified monodromy matrix. Using this as our starting point we extract the first couple of 
the local integrals of motion, and the corresponding time components of the Lax pairs for 
both models. Due to analyticity requirements imposed on the time components of the Lax 
pairs, certain sewing conditions emerge at the defect point. The involution of the charges is 
shown based on the underlying algebra, and the invariance of the sewing conditions under 
the Hamiltonian action is also explicitly checked. 



2 The Landau-Lifshitz model with defect 

We first consider integrable point-like dynamical defects in the context of the isotropic 
Landau-Lifshitz model. We assume periodic boundary conditions and restrict our atten- 
tion to the SU2 classical algebra [23]. Results regarding higher rank algebras may also be 
considered [24]. 

The equations of motion associated with the isotropic Landau-Lifshitz model are given 

as: 

where the vector-valued functions S{x) = {Si{x) , S2{x) , Ss{x)) which describe the physical 
quantities of the model take values on the unit 2-sphere, i.e. S ■ S = 1, and satisfy an SU2 
Poisson structure given by the following Poisson brackets 

Sa{x), Sb{y)} =2teabcS,{x)S{x-y), (2.2) 
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with Eabc being the totally antisymmetric Levi-Civita tensor with value £123 = 1. We note 
that throughout the present section we shall also use the following linear combinations of 
the fields Si{x): 

S^ix) = ^{S^ix)±tS2{x)). (2.3) 
The Hamiltonian and the momentum of the model are given by the following expressions 

n = -- / dx —] , 



4 J_j^ \dx J 

where the prime denotes derivative with respect to x. Using the Hamiltonian, the equations 
of motion can also be expressed as 

'dt 



n, S\. (2.5) 



2.1 The Lax and defect operators 

We may now introduce a single point-like defect of dynamical type in the L-L model at 
the point xq. The starting point for analyzing the system is the derivation of the modified 
monodromy matrix, which is built as a co-action [1] 

r(L,-L,A) = T+(L,xo,A)P(xo,A)T-(xo,-L,A) 

= Pexp(^^^ dx U+(a;)) V{xo, X)Vexp(^ j " dx U-(x)) . (2.6) 

Recall the Lax pair (U, V) for the bulk L-L model 

As indicated in e.g. [1], the defect operator is required to satisfy the same quadratic algebra 
as the monodromy matrix, i.e. 



|P(A)®P(^)}= [r(A-/i),P(A)®P(/i)]. (2. 



In the case of the L-L model, the classical r-matrix is the Yangian solution [25] 

r(A) = -^, (2.9) 
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with V being the permutation operator: V{d ® h) = h ® a. It turns out that the generic 
defect operator satisfying the quadratic Poisson structure with the Yangian r-matrix has the 
following form 

V{\) = \l+(^^^ (2.10) 
due to (2.8) the elements satisfy the 3^2 algebraic relations 

{e, r} = 2^. (2.11) 

Notice that the associated Casimir is given as: C = {^^Y + ^^C,~ hence, the algebra (2.11) 
may be parameterized by two free fields. 

The continuum "bulk" monodromy matrices satisfy the usual differential equation 

^^^^^ = U±(x,A)T±(x,y;A), (2.12) 

and the zero curvature condition is then expressed as: 

i]^{x,t) -Y^\x,t) +\l]^{x,t),Y^{x,t)\ =0, x^xq. (2.13) 

On the defect point in particular the zero curvature condition is formulated as (see e.g. 
[17, 1] for more details) 

^P(xo) = Y+{xo)V{xo) - V{xo)Y~{xo) , (2.14) 

and describes explicitly the jump occurring across the defect point. This will be a major 
consistency check of the prescription followed here. The time components Y^ of the Lax 
pairs for the left and right bulk theories as well as the defect point will be explicitly derived 
subsequently. 

2.2 Local integrals of motion 

We first derive the tower of involutive local integrals of motion. Regarding the left and right 
bulk parts, we consider the usual ansatz for the respective monodromy matrices 

T±(x, y; A) = (1 + W^{x)) e^^^"'^) (1 + W^{y))-^ , (2.15) 

with and being purely off-diagonal and diagonal matrices respectively. We also 
assume that they admit an expansion in terms of the spectral parameter A as 

oo oo 

W^{x, A) = ^ A"l^±(x) , Z^{x, y,X)=Y. ^"^n y) ■ (2-16) 

n=0 ri=— 1 
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Substituting the ansatz (2.15) into the relation (2.12), and sphtting the resulting equation 
into a diagonal and an off-diagonal part one obtains 

W^\]^ - \]^W^ + w^vtw^ - = 0, 



^ = U^ + U^iy±. (2.17) 

The off-diagonal part corresponds to a typical Riccati type differential equation. Solving 
this set of differential equations provides explicit expressions of the W^^Z"^ matrices. As 
for the local integrals of motion, recall that they are provided by the generating functional 

6;(A) = ln(trr(A)) . (2.18) 

Substituting the ansatz (2.15), the latter expression can be formulated as: 

^(A) = lntr[e^"'(^'"«)(l + W+{xq))-^ V{xo) (1 + Vr-(xo))e^"("'"-^»] , (2.19) 

where Schwartz boundary conditions at the endpoints x = ±L have been implemented. 

Solving (2.17) one determines and Z"^ order by order. The first three terms of the 
expansion suffice to compute the Hamiltonian and the momentum of the bulk theories, and 
are found to be 

Oil/X) : W^o± 
C(A°) : -- 



-a\ 1-5, 2S^ 



a J 2S- I + S3 

-a'' 
-a' 



Note that a is obtained by interchanging ^ 5^. Correspondingly, we find that Z^ is 
given by 

"1 00 /* X 

Z±(x,|/,A) = — (x-y)a^ + VA"-^ / dziS+a- + S-a+)W^{z) . (2.21) 

n=i 

At order (9(1 /A) one finds 

Z^, = ^{x-y)a^ . (2.22) 

This term is important, since it provides the leading contribution of as A —i- 0, a fact that 
is used when expanding the generating function Q. The next two orders provide the first two 
physical integrals of motion for the right and left theories. More specifically, one concludes 
that 

{Z^)n = -V^, and (Z±)n = • (2.23) 
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It should be stressed out that the implementation of special boundary conditions at the 
defect point Xq amounts to the emergence of certain boundary type terms in both charges, 
which have to be taken into account. 

Expanding then the generating function Q in powers of A leads to the local integrals of 
motion of the defect L-L model. More precisely, the generating function of the local integrals 
of motion becomes 



GiX) = Z+iX) + Z{1(A) + In [(1 + W+ixo)) ' Vixo) (l + l^-(xo)) 
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(2.24) 



so that the terms Zj^ provide the left and right bulk charges computed above, whereas the 
third term gives the defect contribution. 

Gathering all the information given above we derive the explicit expressions for the 
momentum and Hamiltonian of the system: 



V 



2././" 1 + 53 +2 



dx ^1^2-f^2^1 
XQ 1 + 'S's 



+ I ln(l + Ssixo)) - I ln(l + 53(^0)) - \nV^°\ 



xo) 



n 



where we define 



1 

~4 
1 .X 



r \dx) 



44 \dxJ 



2 ('S'i(xo) - S2{xo)) a+'(a;o) + ^{Si{xo) - S2{xo)) a '1 



Xo 



25^ 



25^ 



(2.25) 



(2.26) 



2^- 1 + 5-3 25- l + S'3 

and are defined analogously. Clearly the variables Si describe the left bulk theory, whereas 
the variables Si describe the right bulk theory. Note also the emergence of some non-trivial 
boundary type terms at the defect point, in both charges. Regarding the defect contributions, 
we have found the following explicit expressions for the first two physical charges 

1 



p(o) 



a-C + + (1 - d+a-)C 
I?« = i [o+(l + a-a+)+ (a+' + (a+)V)e' 

a+a+')a" + C 



with 



Z = 0+(a-r + a^^^ + (1 - a+a-)C 



(2.27) 



(2.28) 



These two terms are the defect contributions in the momentum and Hamiltonian charges 
respectively. For the sake of clarity, we regard from now on a''=,a^ as the fundamental 
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variables, instead of the fields Si, Si. We should also note that these expressions simplify, 
when one considers the sewing conditions on the defect point, as will be transparent in the 
subsequent section. 



2.3 The modified Lax pair 

From the general theory of classical integrable systems, the time components V of the Lax 
pair are given by [1, 23] 

Y{x,X,fi) =t-\X)tTa(raiL,x,X)rabiX,fi)ra{x,-L,X)^ . (2.29) 

In the case where defects are taken into account, one has to compute V in the bulk, as well 
as at the defect point. If the r-matrix of the model is the Yangian, as it happens for the L-L 
model, the corresponding time components are simplified and for a single point-like defect 
are expressed as [1] 

V+(x,A,/i) = T+{x,xo)V{xo)T-{xo,-L)T+{L,x) 

X — fjL 

_ t-1 _ 

V (x, A,yu) = - T {x,-L)T^{L,xq)'D{xq)T [xq^x) 

X — jj, 

V+(xo,A,/i) = P(a:o)T-(a:o,-L)T+(L,Xo) 

X — jj, 

V-(xo,A,/i) = T-(xo,-L)T+(L,xo)I^(xo) . (2.30) 

A — /i 

We provide explicit expressions for the first three orders of the V operators below. It should 
be clear that the following expressions hold for the left bulk part 



l + S _ 1/1 + 53 2S' 
~2ir ~ ~2il [ 2S+ 1 - ^3 



0{X'): V^^^(x,/i) = -^ = --( ^ ^1 (2.31) 

fi\ S3S+' - S'^S+ S~S+' - S+S-' ) ' /i 



Similar expressions hold for V^^, but with S** — )■ 5**. At order 0{X) one naturally finds the 
time component associated with the Hamiltonian. For the sake of clarity, as noted before, 
it is more useful to write down the operators in terms of the fields a^, defined in (2.20) 

0(A°): %^{x,f^) = -y^^(\ t) (2.32) 
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At the defect point, we have the following expressions at the first order 



Vj))(xo,/i) = 

V,-„,K.,) = -~{ J . (2.33) 



Gluing the bulk and defect V operators at the defect point leads to sewing conditions 
that relate the bulk fields at the defect point with the defect fields. In particular, setting 

V+)(x->xo)=V+)(xo) , (2.34) 

leads to a unique sewing condition for the four matrix entries 

C(°) : ^+ - fl-a+r - {a' + a"^)^ = . (2.35) 

Solving this condition with respect to and substituting the solution back to the defect 
contribution of the first charge (2.27) simplifies the latter one significantly, which can be 
written in a very compact form as 

p(o) = a-^- + r • (2.36) 

The sewing condition emerging when gluing the V operators at the defect point from the 
left is related with eq. (2.35) by complex conjugation. 

At the next order, the modified V operator at the defect point from the right is 



with 



Qu = (^a^a+' + (a"' - a^)a+j + (^a"a+ - a^a+' - 2^ CC 

+ (a+' - a+) + {a^a+ - 1 + a'"a+ - a~a'+) i^)^ - (a^C? (2.37a) 

+ (^a+(-l + a-a+-l + a+a-') (r)^-(a^ + a^')(a"D^ (2.37b) 

= {a- fee + (2a+a- - 2a-a+')^+e - (a" + 0"')^+^ 

+ (a- - a+{a~y + {a~ya+' - a~') (e)' + («+' - a+){t? (2.37c) 
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= -{a-a+' + a+{a- + a-'))er + {a~)\a+ + a+')rC-{a^tf 

(2.37d) 

Equating the first order Lax operators at the defect point, V^-)(x — )■ Xq) = V^^(xo), yields 
a unique second sewing condition, C^^\ for all four entries of the matrices, which is rather 
complicated and long. However, using the previously found sewing condition of the lower 
order, C^^\ the new condition simplifies greatly and assumes the form 

a+(r)+)^(a-r + D + a"a"^'r + + a"'r) + (a"' + = , 

or, in other words, 

C(^) : (a- - a+) + (a-a+)' T + (a^ + a^)' T = , (2.38) 

since the factor a~^{d~^y in front can not be zero and the next term inside the parenthesis is 
just T)^'^\ We observe that the sewing condition at this order involves not only the left and 
right fields, as well as the defect fields, but also the derivatives of the bulk fields (see also 
[1, 2]). It has been explicitly checked that the extracted sewing conditions are as expected 
compatible with the Hamiltonian hierarchy, confirming the relevant generic statement proven 
in [1] regarding the invariance of all sewing conditions under that Hamiltonian action. Note 
also that the commutativity of the momentum and Hamiltonian is valid as expected, but 
on-shell, i.e. the sewing conditions must be taken into account, as it happens in the case of 
the sine-Gordon model with defect (see also [2]), 

n, P| = 0. (2.39) 

It is worth noting that taking into account the second sewing condition (2.38), the first 
order defect contribution simplifies significantly; it is eventually given by the following simple 
expression 

= + . (2.40) 

a+ — a 

It is interesting that on-shell the defect fields C,^ cancel out completely from the defect 
Hamiltonian. 

Finally, the equations of motion for the left and right bulk theories are given by (2.5), 
for the right theory the variables are simply Si. The time evolution on the defect point is 
expressed as: 



r = I (a - I)(i)a+r)+)r + ^ (V^'^^+il - a-a+) - c)r 

e+ = ^ ( - B + I?«a-0+)r + 4 ( - V^'^-O^il - a-a+) + 0)^"^ 
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t = (a - p(i)a+0+)e^ ^^'^^"^^)^" ' (2.41) 

where we define 

A = a+' + (a+)^a+', B = (a+ - a+')a"a+' + t)+a"', 

C = a+'(a+ + a~) + a+(a+a+'a" - a+' - 0+a"') . (2.42) 

The equations above, relevant to the defect point, emerge either from the Hamiltonian 
description through 

C = {h, C] , (2.43) 
or the zero curvature condition on the defect point (2.14). 



3 The Faddeev-Reshetikhin model with defect 

The F-R model was introduced in [19], as a convenient modification of the familiar non- 
ultra local PCM model. Despite the different underlying algebras, it was shown that the 
equations of motion of these two models coincide, hence they may be considered as equivalent. 
Moreover, the F-R model is easier to quantize, due to the absence of non-ultra local terms 
in the algebra. 



3.1 The model 

The Lax operator of the model is the same with the familiar PCM one and has the form 



,x A Jo(x) + Ji(x) 
U(x,A) = ' ^^-^^ 



with Jj being currents taking values on the Lie algebra g, corresponding to a Lie group G. 
In the present work, we restrict ourselves to the simple case of the SU (2) group. We choose 
the Pauli matrices (a"'', a", \cr^) as the generators of the SU2 Lie algebra. We also define the 
currents 

5 = J^ = ^, T = j- = :^, (3.2) 

so that the Lax operator may also be written as 

Expanding in the Lie algebra generators, the currents are expressed in the following compo- 
nent form 

S = S^t" = S'^ + S^a+ + S-a- , (3.4) 
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and similarly for T. The time component of the Lax pair is given by: 

V(.,A) = -^ + ^. (3.5) 

The Lax operator satisfies the fundamental Poisson bracket relation 

V{x,\)®V{y,i^i)] = [r(A-/i), U(x, A) ® I + I ® U(x, /i)! 5(x - , (3.6) 



which reflects the absence of non-ultra-local terms in the level of the algebra for S and T. 
In particular, the latter ones satisfy the following algebraic relations 



[S\x), T\y)] = 

[T^{x),T\y)} = r\T\x)5{x-y). (3.7) 

As mentioned above, we restrict our attention in the simplest case, i.e. the SU (2) group. The 
underlying algebra consists of two independent copies of the SU2 Lie algebra and the structure 
constants are proportional to the usual Levi-Civita tensor {f^\ = 2ieabc)- Moreover, S and 
T may be interpreted as three- vectors and their lengths can be fixed as: \S\ = \T\ = C = 
constant] i.e. they are simply the SU2 Casimir operators. Note that the latter restriction 
(C = constant) is known as the Virasoro constraint in the string theory frame (see e.g. [20]), 
rendering essentially the familiar PCM model ultra-local. 



3.2 Local integrals of motion 

Let us first briefly review the bulk theory. From the Lax operator of the F-R model, it is 
manifest that there are two expansions that one may obtain, around two singular points. 
The local conserved charges, such as the momentum and the Hamiltonian are then linear 
combinations of the quantities derived. We present details for the A = l, (y = A — 1 = 0) 
expansion, the results for the A = — 1 are similar and are omitted. In particular, for the bulk 
parts we have 

oo 
n=0 

and 

oo 
n=-l 

For the first order we have found for W 

«o = ^{S' T C), -ao = bo = ^("^^ ± C) , (3.10) 
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with 

C = >y{S-y + A{S-){S+) = 2tT{S^) = 2tT{T^) = ^(T-)2 + 4(T-)(T+) . (3.11) 

Setting the Casimir C equal to one amounts to coinciding with the results of the L-L model. 
In the next order 

± ai (2C) = -T- - ao + al T+ + 2a[) 

±h{2C) = T+ -boT' -blT- -2b'Q . (3.12) 

Observe also that 

bi = —ai + I , (3.13) 

so that everything can be written in terms of and Oi, which may now be regarded as the 
fundamental fields. We carry on the computations for the Z elements. The results we have 
found are the following: 

2 = +S' + 2bo S' = ±C 
2 6'_^ = -S' + 2ao S+ = tC 
4 7^, = +T' + 2bo + 46i S~ 

Ad'o = -T' + 2ao T+ + Aai S+ . (3.14) 
It is manifest that the first integrals are trivial since the corresponding densities are constant 

X_i=X_i = LC. (3.15) 

Henceforth we choose to consider C = 1, as well as the first of the two signs appearing 
in all the expressions above. Let us first derive the expression for the bulk theory with 
Schwartz type boundary conditions. Substituting the exact expressions for a^'s and 6j's we 
obtain the charges in terms of the fields. Bear in mind that some boundary terms vanish due 
to Schwartz boundary conditions, that will have to be taken into account when considering 
non-trivial boundary conditions at the defect point. Eventually, one is left with 



1„ 1 



L 



Vs-- dxti{ST) 



^Vt + - I dxtT{S T) , (3.16) 



and 



2 2 

, SiS'-S[S2 _ , T{T'-T[T2 , 

Vs = i \ dx \' \ VT = i dx \^ ' ^ . 3.17 

Taking their sum and their difference yield the Hamiltonian and momentum densities of the 
system respectively. More specifically, 

'P = \Ps + \Vt 
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n = ]^Vs-\Vt- j ^dxii{ST) . (3.18) 



The above results hold for the right and left bulk parts of the F-R model. Recalling the 
generic expression for the generating function (2.24), we find the defect contributions at all 
orders. 

Now consider the system in the presence of a point-like integrable defect at the position 
Xo, described by the generic L- matrix (2.10). The expressions of the momentum and the 
Hamiltonian emerging from the expansion of A around the values ±1, are given as: 

- iln((l + 53)(l + ^3))(xo)-ln2J(xo), (3.19) 

- lln((l+T3)(l + f3))(xo)-ln*ij(xo), (3.20) 
where we define 

^ = r(i + 4^0) -ao+e+ + &or- 1 + 4^0, (3.21) 

where Oq , 69 defined as in (3.10) with the S fields, , 6q are defined with the S fields, 
and 

Ts-1 ^_ 1 - T3 . , f3 - 1 ^ , 1 - f3 

an = — Or. = — — , On = — = — , On = — = — . (3.22) 

2T~^ 2T~ 2T+ 2T~ 

The fields S*, T correspond to the left theory whereas the fields S, T correspond to the right 
theory. The complete momentum and Hamiltonian of the system are given by the following 
expressions: 

V = Xo+Xo 

n = Xo-Xq. (3.23) 



3.3 The modified Lax pair 

The next step is the derivation of the associated time components of the Lax pairs for the 
bulk theories and at the defect point, from left and right. Requiring continuity around the 
defect point the relevant sewing conditions are extracted. The bulk quantity associated to 
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the Hamiltonian is given by (3.5). The V operator at the defect point from the right, from 
the expansion around A = 1 has the following expression: 



1 



e + 1 + 



(A - 1)5J + b,i-e + 1) -ate - 4boi-e + 1), 
while from the expansion around A = — 1 one obtains 



(3.24a) 



1 



«o ^ r 



(A + 1)^1 \e + &o (-e - 1) -ate - atboi-e - K 
where 2J, QJ are defined in (3.21). Similar expressions hold for V~, in particular: 



(3.24b) 



and 



e + 1 - 



atie 



(A - 1)23 \-atb,e + boie + 1) &o r - (-e + 1: 



1 



4r 



a+(e + i) + r 



(3.24c) 



(3.24d) 



(A + 1)2J \-atb^e + feo (e - 1) &o r + a+b^ie + 1) 
Due to the continuity requirements, the corresponding sewing conditions arise, from the 
expansion around A = 1: 



'0 ~ bo 



0, 



c^: e-btbor-{bt + b,)e-m 

and from the expansion around A = — 1: 

: e - btb^r - {bt + bo)C - {bt - &o ) = • 



(3.25) 



(3.26) 



Commutativity between the momentum and the Hamiltonian of the F-R system is shown 
via the algebraic relations of the defect degrees of freedom, and as expected: 



{n, v] 



0. 



(3.27) 



Furthermore, as in the case of the L-L model, the compatibility of the sewing constraints 
with the charges in involution is explicitly checked and confirmed. 

Finally, by means of the underlying Poisson structure, or the zero curvature condition, 
one may derive the equations of motion for the bulk theories 



S{x,t) = -S'{x,t) - S{x,t), T{x,t) 
f{x,t) = T'{x,t)- \T{x,t), S{x,t) 



x^ Xo 



(3.28) 
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Similar equations of motion hold for the fields of the right theory S, T, which are omitted 
here for brevity. The time evolution of the defect degrees of freedom reads as: 



/a a \ , / c c 
-2( ^1 



where we define 



a = — ag , b = 6q , c = 1 + a^b^ 

a = — dg , b = feg", 6 = 1 + 0060'. (3.30) 

Note that expressions (3.29) arise either from the Hamiltonian or from the zero curvature 
condition at the defect point. With this, we complete our analysis on sigma models in the 
presence of dynamical point-like defect. In the last section we summarize our results and 
discuss some interesting extensions that will be analyzed in future works. 



4 Discussion 

We have applied the methodology developed in [1] for the Landau-Lifshitz and the Faddeev- 
Reshetikhin models. Due to the existence of the classical r-matrix, we have been able 
to formulate a modified monodromy matrix, and hence identify the hierarchy of Poisson 
commuting Hamiltonian alongside the associated Lax pairs. Certain sewing conditions com- 
patible with the Hamiltonian action have also been identified. 

In the present and previous studies dynamical type defects have been considered, the so 
called type-II defect according to the terminology of [6, 8]. This type of defects is associated 
to a Poisson structure strictly related to the point-like defect, hence time evolution of the 
defect degrees of freedom occurs. The next natural step would be the investigation of the 
type-I or non-dynamical defects. Study of non-dynamical defects may lead to Backlund 
transformations associated to the defect point as pointed out in [6, 8]. It would be quite 
interesting to check whether these Backlund transformations could naturally emerge through 
the proposed algebraic formulation. 

We have restricted our attention so far to point-like defects, although in principle spatially 
extended impurities may also be investigated within the proposed scheme. Moreover, at the 
quantum level one may study the associated Bethe ansatz equations, and hence extract 
significant information regarding the interaction of particle like excitations displayed by 
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such models with the defect. All the above are highly intricate issues and hopefully will be 
addressed in forthcoming publications 
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